BALLISTIC  RESEARCH  LABORATORIES  REPORT  NO .  693 


TABLES  FOR  THE  PIDDUCK-KENT  SPECIAL  SOLUTION 
FOR  THE  MOTION  OF  THE  POWDER  GAS  IN  A  GUN 


by 

John  P.  Vinti 
and 

Sidney  Kravitz 


ABERDEEN  PROVING  GROUND,  MARYLAND 
JANUARY  1949 


Form  SF298  Citation  Data 


Report  Date  D  r 

("DD  MONYYYY")  eport  ype 

01011949  N/A 

Dates  Covered  (from...  to) 

("DD  MONYYYY") 

Title  and  Subtitle 

Tables  for  the  Pidduck-Kent  Special  Solution  for  the  Motion  of 
the  Powder  Gas  in  a  Gun 

Contract  or  Grant  Number 

Program  Element  Number 

Authors 

Project  Number 

Task  Number 

Work  Unit  Number 

Performing  Organization  Name(s)  and  Address(es) 

Ballistic  Research  Laboratories  Aberdeen  Proving  Ground,  MD 

Performing  Organization 

Number(s) 

Sponsoring/Monitoring  Agency  Name(s)  and  Address(es) 

Monitoring  Agency  Acronym 

Monitoring  Agency  Report 

Number(s) 

Distribution/Availability  Statement 

Approved  for  public  release,  distribution  unlimited 

Supplementary  Notes 

Abstract 

Subject  Terms 

Document  Classification 

unclassified 

Classification  of  SF298 

unclassified 

Classification  of  Abstract 

unclassified 

Limitation  of  Abstract 

unlimited 

Number  of  Pages 

50 

AMSRL-CS-IO-SC  (380) 


28  July  2000 


MEMORANDUM  FOR  Chief,  ARL  Technical  Library,  ATTN:  Ms.  L. 
LeTendre 

SUBJECT:  Distribution  Statement  for  BRL  Report  No.  693 


1.  References: 

a.  Memorandum,  AMSRL-CI-LP,  undated.  Subject:  Review  for 
Change  in  Classification  and  Distribution,  ARL  Report  No.,  copy 
enclosed. 

b.  Ballistic  Research  Laboratories  Report  No.  693,  "Tables 
for  the  Pidduck-Kent  Special  Solutions  for  the  Motion  of  the 
Powder  Gas  in  a  Gun",  by  J.  P.  Vinti  and  S.  Kravitz,  January  1949, 
copy  enclosed. 

2 .  The  correct  distribution  statement  for  the  referenced  report 
is : 


Approved  for  public  release;  distribution  is  unlimited. 

3 .  Request  that  you  mark  all  of  your  copies  of  the  report  with 
this  distribution  statement.  We  have  notified  the  Defense 
Technical  Information  Center  about  the  distribution  statement  for 


TABLE  OF  CONTENTS 

Page 

Abstract  1 
Section  I  Elementary  Corrections  for  the  Motion  of  the  Powder  Gas  3 
Section  II  The  General  Problem  of  the  Motion  of  the  Powder  3 
Section  III  The  Lagrange  Problem  4 
Section  IV  The  Pidduck-Kent  Special  Solution  4 
Section  V  An  Improved  Derivation  of  the  Pidduck-Kent  Special  Solution  5 
Section  VI  Discussion  of  the  Above  Results  13 
Section  VII  The  Methods  of  Calculation  and  Tabulation  15 
Section  VIII  Discussion  of  the  Tables  and  Their  Accuracy  22 

Appendix  A  Tables  for  the  Pidduck-Kent  Solution  27 
Appendix  B  Graphical  Representation  30 
Appendix  C  The  Evaluation  of  Certain  Integrals  34 
Table  of  Symbols  ^6 


«■  1 


Atrstract 


A  discussion  is  first  given  of  the  general  problem  of  the  flow 
of  gas  in  a  gun,  of  the  simplified  problem  known  as  the  Lagrange 
problem,  and  of  the  Pidduck-Kent  special  solution  of  the  differential 
eauation  of  the  Lagrange  problem.  This  special  solution  is  commonly 
supposed  to  be  the  limiting  solution  of  the  latter.  An  improved 
derivation  is  then  given  of  the  Pidduck-Kent  solution,  the  improve¬ 
ment  consisting  in  accounting  for  gas  imperfection  by  means  of  a 
co volume  correction. 


The  Pidduck-Kent  solution  is  then  put  into  a  form  suitable  for 
easy  calculation.  To  summarize  the  results,  let  the  following  sym¬ 
bols  denote  certain  quantities  as  calculated  from  the  Pidduck-Kent 
solution,  viz 

pQ  =  breech  pressure 

pb  -  pressure  at  the  base  of  the  projectile 

P  -  space-mean  pressure 

Wp  =  kinetic  energy  of  the  projectile 

W_  -  kinetic  energy  of  the  powder 
c 

Then,  according  to  the  Pidduck-Kent  solution, 

-jr-  =  •  +€/s 

■/a,  -  + «5ia 

f  ^  c/ 

where  6  denotes  the  charge-projectile  mass  ratio  An,  H  the  polytropic 
index  •/(If- 1),  (where  if  is  the  effective  ratio  of  specific  heats  that 
takes  heat  loss  into  account) ,  and  a  a  parameter  characteristic  of  the 
Pidduck-Kent  solution.  * 


The  quantity  &  is  given  by 

i-ri«[,+aCB( 


1  +  c^wVj 
i+c,n  /J 


The  quantity  o< ,  depending  only  on  6  ,  may  be  found  by  linear  inter¬ 
polation  in  Table  I  of  Appendix  A.  The  quantity  ,  also  depending 
onJ.y  on  €  ,  may  be  found  by  linear  interpolation  ^n  Table  II  of 
Appendix  A.  The  quantity  C.,  depending  weakly  on  both  £  and  H  ,  may 
be  found  by  two-way  linear  Interpolation  in  Table  III  of  Appendix  A 
or  by  reading  from  the  graphs  or  the  contour  map  in  Appendix  13 .  Ordinarily 
reading  from  the  graphs  is  the  easiest  and  most  accurate  of  the  three 
methods  of  finding  c,.  With  these  aids  one  may  expect  to  obtain  Pidduck- 
Kent  values  (as  distinguished  of  course  from  true  values)  with  the  follow¬ 
ing  maximum  errors  for  the  range  £  =  0  to  10,  if  -  1.2  to  3: 
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I  PART  IN  5000 
I  PART  IN  6700 


I  PART  IN  5000 


Section  I.  Elementary  Corrections  for  the  Motion  of  the  Powder  Gas 


In  cases  where  the  initial  mass  £  of  the  powder 

is  small  compared  to  the  mass  m  of  the  projectile,  it  is  customary 
to  correct  for  the  motion  of  the  powder  gas  by  means  of  certain 
very  simple  formulas.  These  are; 

R-P.O+%)  (l) 

P-jpi+«/3)  <2) 

wc“  (VsHIfemV*)  (3) 

where  c/m  (4) 

and  pQ  denotes  the  pressure  at  the  breech,  pfc  the  pressure  at  the 

base  ,  of  the  projectile,  p  the  spaee.T-mean . pressure  behind  the  pro¬ 
jectile,  V  the  velocity  of  the  projectile,  and  V  the  kinetic 

energy  of  the  powder  plus  powder  gas.  These  formulas  are  derived 
on  the  assumptions  that  the  cross-sectional. area  is  uniform  all 
the  way  from  the  breech  to  the  base,  of  ...the  projectile  and  that  the 
powder  grains  and  powder  gas  together  form  a  fluid  of  uniform  density. 
Under... these.,  conditions  the  velocity  of  this  fluid  varies  linearly 
ffoifl  the  value  zero  at  the  breech  to  the  value  V  at  the  base  of  the 
projectile.  The  correction  for  the  kinetic  energy  of  the  charge, 
vifc  addition  of  one-third  the  mass  of  the  charge  to  that  of  the  pro¬ 
jectile,  is  then  similar  to  that  for  the  kinetic  energy  of  a  spring 
(without  waves)  to  which  a  oob  is  attached. 


Section  II.  The  General  Problem  of  the  Motion  of  the  Powder. 

Actually  the  problem  is  somewhat  more  complicated  than  the 
simple  theory  would  indicate.  Gas  friction  is  not  taken  into 
adeoUht  at  all  in  the  above  formulas  for  the  pressure  ratios. 
Chambrage,  the  drop  in  cross-sectional  area  on  going  from  the  ’ 
powder  chamber  to  the  bore,  introduces  a  non-uniformity  into  the 
cross-sectional  area,  a  non-uniformity  which  will  probably  become 
much  more  serious  as  one  goes  in  the  direction  of  guns  of  much 
higher  velocity.  The  powder  gas  is  compressible,  so  that  there 
actually  must  be  a  variation  of  density  with  position  as  well  as 
time.  It  is  not  really  correct  to  lump  the  powder  grains  and  the 
powder  gas  together  as  a  single  fluid.  Finally,  the  burning  of  the 
powder  furnishes  energy  to  the  gas,  so  that  the  expansion  is  not 
adiabatic. 
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Section  III.  The  Lagrange  Problem. 


When  Lagrange  set  himself  the  problem  of  solving  for  the 
motion  of  the  powder  gas  in  a  gun,  he  was  able  to  formulate  it 
in  a  relatively  simple  fashion,  because  the  propellant  of  his 
day  was  simply  loose  black  powder.  In  such  a  case  it  is  not 
too  bad  to  assume . instantaneous  combustion.  He  therefore  assumed 
that  he  had  a  cylinder,,  closed  at  one  end  by  the  breech  arid  at 
the  other,, variable,,  end  by  the  base  of  the  projectile  and  filled 
initially  with  hot  gas  of  uniform  pressure,  density,  and  tempera¬ 
ture.  .His  problem  was  then  to  calculate  the  subsequent  states 
of  the  system,  the  projectile  being  initially  at  rest. 

Lagrange's  problem  has  received  its  most  complete  treatment 
by  Love  and.  Pidduck*,  who  applied  their  results  to  the  calcula¬ 
tion  of  p  /ph  and  V  /V  versus  travel  for  the  case  of  a  150  mm. 
gun.  (W  °aenotes  the  ^kinetic  energy  of  the  projectile.)  Pidduck 
found  that. these  ratios  oscillated  but  seemed  to  approach  certain 
limiting  values,  whioh.  values  corresponded  to  a  certain  special 
solution  of  the  differential  equation  of  motion  of  the  gas.  The 
rigorous,  solution  is  of  a  wave  character,  involving- rarefaction 
waves  travelling  . back  and  forth  between  breech  and  projectile, 
and  satisfies, the... initial  conditions  of  the  Lagrange  problem.  The 
special  solution,  on  the  other  hand,  is  of  a  non-wave  character 
and  does  not  satisfy  the  initial  conditions  of  the  Lagrange  pro¬ 
blem,  but  corresponds  to  an  initial  non-uniform  distribution  of 
pressure  and  density.  From  observation  of  the  computational 
results  for  PQ/P^  Pidduck  and  all  later  investigators  have  sus¬ 
pected  but  not  proved  that  the  accurate  solution  approaches  the 
special  solution  in  the  limit  of  large  travel. 


Solution  IV.  The  Pldduck-Kent  Special  Solution. 


The  special  solution  has  also  been  derived  by  Kent**  and 
applied  by  Hirschf elder  et  al***  to  the  computation  of  the  ratios 
and"W m  the  limiting  case  of  sufficiently  small 

values  of  £  ,  the  Pidduck-Kent  special  solution  gives  values  for 
the  above  ratios  that  agree  with  the  uniform-density  values  of 
equations  (2)  and  (3) .  The  deviations  from  the  uniform-density 
values  become  important  in  the  case  of  large  values  of  £  ,  i.e. 
for  high-velocity  guns.  As  we  have  said  above,  the  effect  of  gas 


*Love  and  Pidduck,  Phil.  Trans.  Roy.  Soc.  222,  167  (1922) 

**R . H .  Kent,  Physics  7,  319  (1936) 

***Hirschf elder,  Kershner,  and  Curtiss,  NDRC  Report  A-142 
Hirschf elder,  Kershner,  and  Sherman,  NDRC  Report  A-204 
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friction  also  becomes  more  important  for  high  velocities,  as  does 
the  effect  of  chambrage.  The  latter  becomes  more  important  because 
the  use  of  a  high c/m ratio  necessitates  the  use  of  a  powder  chamber 
much  fatter  than  the  bore.  For  small  values  of  e,  for  which  the 
assumption  of  uniform  density  may  be  made,  the  effect  of  chambrage 
may  be  estimated  by  the  methods  of  BRL  Reports  No.  and  No.  3* **. 

For  large  values  of  e  it  is  doubtful  that  those  methods  would  apply. 

Although  the  solution  of  the  hydrodynamical  problem  thus  in¬ 
volves  a  number  of  complications  when  we  go  to  large  values  of  e, 
it  appears  worth  while  to  isolate  the  effect  of  non-uniformity  of 
gas  density  by  publishing  adequate  tables  of  the  Pidduck-Kent  special 
solution.  These  tables  may  not  eventually  be  useful  for  practical 
calculations,  but  they  should  be  of  considerable  help  in  the  analysis 
of  experimental  results  on  high-velocity  guns,  if  only  to  tell  us,  by 
comparison  of  measured  pressure  distributions  with  the  Pidduck-Kent 
solution,  how  important  are  the  neglected  factors. 


Section  V.  An  Improved  Derivation  of  the  Pidduck-Kent  Special  Solution 

In  the  Lagrange  problem,  of  which  this  is  a  special  solution  of 
the  corresponding  differential  equation***,  it  is  assumed  that  the 
chamber  and  bore  form  a  cylinder  of  uniform  cross-section,  the  powder 
being  all  burned  initially.  We  follow  the  notation  of  Kent  in  the 
above-mpntioned  paper. 

Thus 

x  =  initial  distance  of  any  gas  particle  from  the  breech 
y  =  distance  at  time  t  of  the  same  gas  particle  from  the  breech 
p  a  density  of  the  gas 
p  =  pressure  of  the  gas 


The  subscript  zero  denotes  initial  values.  A  subscript,  on  a 
partial  derivative  denotes  the  variable  which  is  kept  constant  in 
the  differentiation.  The  equation  of  continuity  is  then 


(5) 


*J.P.  Vinti,  BRL  Report  No.  309,.  ’’The  Equations  of  interior 
Ballistics",  Chap.  IV. 

**N.F.  Ramsey,  Jr.  (with  J.R.  Lane),  BRL  Report  No.  3,  "Analysis  of 
Pressure-Time  Curves  for  the  Gerlich  Rifle" 

***It  is  not  a  solution  of  the  actual  Lagrange  problem,  because  the 
initial  conditions  are  not  satisfied. 
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The  equation  of  the  state  of  the  gas  is 

r(jH)  -r.t 


(%) 

Cl) 


■where  n  denotes  the  specific  covolume  of  the  powder  gas,  T  the 
absolute  temperature,  and  R^  the  gas  constant  per  unit  mass.  If 


denotes  the  ratio  of  specific  heats  of  the  gas,  the  equation  of 
an  adiabatic  is  . 


?(; 


O) 


Kent  further  assumes  that  "the  original  pressure  and  density  dis¬ 
tribution  follows  adiabatieally  from  an  antecedent  regime  of 
uniform  pressure  and  density."  Under  these  circumstances  (9) 
becomes 


—  CONSTANT  k 


The  derivation  follows  through  equally  well  if  we  replace  the 
true  y'  by  a  larger  effective  value  V  for  the  ratio  of  specific 
heats,  a  device  that  affords  a  means  of  taking  heat  loss  into 
account.*  Furthermore  it  is  not  necessary  to  take  the  right 
side  of  (10)  to  be  a  constant;  it  suffices  for  the  following 
derivation  to  take  it  as(k)(lir(t)),**  where  y(t)  is  a  function 
of  time  alone.  Kent  suggests  that  by  such  a  means  one  can 
allow  partially  for  the  burning  of  the  powder.  It  is  doubtful, 
however,  that  any  such  simple  artifice  will  be  of  use  in 
attempting  to  treat  the  case  where  the  powder  is  still  burning. 
Indeed  the  continuity  equation  (5)  and  the  dynamical  equation 
(6)  do  not  hold  in  that  case,  unless  the  density  p  be  taken  as 
the  density  of  a  "fluid"  composed  of  powder  gas  and  solid  powder 
together,  in  which  case  the  above  relations  between  pressure  and 
density  do  not  hold. 


*  See  J.  P.  Vinti,  BRL  Report  No.  307,  "The  Equations  of 
Interior  Ballistics" 

J.  P.  Vinti,  BRL  Report  No.  402,  "Project  for  a  New  Table 
for  Interior  Ballistics  for  Multiperforated  Powder." 

J.  P.  Vinti  and  Jack  Chemick,  BRL  Report  No.  625,  "Interior 
Ballistics  for  Powder  of  Constant  Burning  Surface." 
♦♦Suggested  by  Kent,  loc.  cit. 
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The  derivation  of  the  special  solution  for  an  imperfect  gas 
has  been  given  partially  by  Kent,  who  showed  that  the  solution  is 
the  same  as  that  for  a  perfect  gas  through  powers  of  e  =  c/m  as 
high  as  the  third.  In  this  section  we  show  that  this  restriction 
can  he  removed,  the  formulas  obtained  being  entirely  independent 
of  the  value  of  l|  .  Pidduck*  has  outlined  a  proof,  but  we  believe 


it  desirable  to  put  the  complete  derivation  on  record. 

Following  Kent,  we  introduce  the  variable 

I7'tJ0dx 

and  attempt  to  find  a  special  solution  of  (6)  in  the  form 

(1-2O 

/Vt  t  =  o  we  have  y  =  x  so  that 

fcx)  (f)(0)  =  X~f].CV*C*Wx  ^ 

From  (12)  it  is  clear  that  either  f  or  <|>  may  contain  an  arbitrary 
factor,  so  that  we  are  free  to  choose  the  function  <}>  in  such  a  way 
that  <Ko)  =  1.  Then 

£Cx)  =  x-  rjf//’e>cx)dx  oh) 

Froia  C^)j  l"),  fWb  (iz)  WE  FIND 

f'tX)$(t)  OS) 


Introducing  the  new  variable 

*  F.  B.  Pidduck,  Journal  of  Applied  Physics  144  (1937) 
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■we  have 


A  -=  web 


where  f »  is  short  for  f ' (x) 


Now 


so  that 


f  (x)  -  i  -  r)/0(x) *]), 

=w»  =  W(t> 


cn) 


CIS) 


Making  use  of  the  fact  that 


WE  HAVE 


05) 


Using  the  continuity  equation  (5)>  we  then  find 


On  writing 


(22) 


In  tern\S  of  w,  (io)  becomes 

p  =  Kv/'l|r(t)  -  Kw0y<|)  *{jr( t) 

From  Ofc),  (2.1),  and  (11) 

w0f  f  (t)  =  -KyWo'"’  (jf'y  ^  (23) 

Separating  variables,  we  obtain 

—  J  ^  . —  ^  WQ  J3  A  CONSTANT  fcVb 

'  '  “  f  df 


Integration  of  (24)j  with  use  of  the  conditions  f 


AND  W0 

-  wo(o)  at  x®  o  ,  g  iv  es 

y-i  .  .  y-*,  v  B(^-0  f Z 

w©  "  w0  (o) 

(25) 

OR 

i 

w„(£)  3  w„(o)(i-af  * 

(16) 

WHERE 

asB(V-')  w0‘  *(o) 

ZK3 

(27) 
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Let  xb  denote  the  initial  distance  from  the  breech  to  the  base 
of  the  projectile,  A  the  uniform  cross-sectional  area,  and  abbre¬ 
viate  f(xb)  by  f.  .  Then,  using  the  fact  that  when  the  powder  is 
all  burned  the  total  mass  of  the  gas  is  equal  to  the  original 
powder  mass  c,  we  have 


j£=  A(fbf d£  =  W.J£ 

If  we  now  use  (26)  and  (28),  introduce  the  dunmy  variable 
defined  by  f  =■  fJL  ffc,  and  use  the  abbreviations  ' 


(2 .8) 


a, 


we  find 


c-  =  Aw ,co)£.S 


Y 


M) 

(30) 

(3|) 


At  the  base  of  the  projectile  we  have  from  (12)  and  (22) 

(ff't)  (U) 

fk=  Kw/(£k)(f'sT|r(-t)  (33) 

If  m  denotes  the  mass  of  the  projectile  (increased  perhaps 
by  some  factor  a  little  greater  than  unity  if  one  wants  to  allow 
for  bore  friction),  the  equation  of  motion  of  the  projectile 


then  becomes 


-V 


Then 


£  0"  =  Ark  =  M  w/(£b)<t>"  TIT 

<t>"(p  _  AKw/(£k)  _  o 

If  m£b 

where  the  constant  B  is  the  same  as  in  (24)  and  (27).  Comparing 
(35)  with  (27),  with  the  use  of  (26)  and  (29),  we  find 


(3t) 

(35) 


C3fc) 

(37) 


Dividing  (36)  by  (31),  we  then  find 

g  =  i^  a.  (i-a.)''r/<ir'l)S 

^-1 

If  ^  and  e  are  given,  the  parameter  aQ  may  be  calculated  from  (37), 
The  Ratio  ?°/l. 


From  ( 22 ),  (26),  and  (29)  the  ratio  of  breech  pressure  p  to 
base  pressure  pK  is  given  by  o 

cz&) 


% 
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The  Kinetic  Energy  Wp  of  the  Powder -Gas; 
Since  the  gas  velocity  ( 


x ,  we  have 


pVTTINO  dx-df/f'^  F30M  (|£);  ~ po 

^*1 £  '  =  W0  PROM  Cl  8 ^  WE  OBTAIN 


<3j) 
FROM  (5*), 


Wf=|rtf,,vvcf2df  .<*•> 

At  the  base  of  the  projectile  the  gas  velocity  equals 

the  projectile  velocity  V.  Use  this  fact,  insert  (26)  into  (40), 
replace  f  by  u.£b  ,  use  (29),  and  replace  AweCo)fb  by  c/$  accoro!H6 
TO  (31).  One  fiHOS 


The  Ratio  r/fb _ 


We  next  investigate  the  ratio  of  the  mean  pressure  p  to 
the  base  pressure  pb  ,  where  p  denotes  the  mean  pressure  that 
should  be  used  in  the  equation  of  state.  In  interior  ballistics 
the  equation  of  state  is  used  to  eliminate  the  temperature,  which 
occurs  in  the  expression  for  the  internal  energy  of  the  gas.  To 
write  down  the  total  internal  energy  we  need  to  find  a  mass  average 
for  the  specific  internal  energy  cvT  .  Now 

fvT  '•))  (42) 

The  total  internal  energy  Et  is  then  given  by 


E‘  =  *£' ?  (7  d* 

The  value  used  for  p  in  interior 
average  value 


The  appropriate  mean  pressure  p 
is  then  given  by 


(43) 

ballistic  calculations  is  the 

(44) 

for  use  in  the  equation  of  state 


(45) 


NOW  FROM  (»)  AND  (12)  WE  HADE 


M  ■»b-lJeV‘>‘*)dx  ”  V1c/a 


(46) 
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Then 


(48) 


(♦9) 


From  (4S)  ano  (4*7)  -»  , 

Now  POT  =  /°J  *^X  “  ^  *  f*"  ”"  W*  ' 

-L.-.rjsw'1,  ^|(0  use  (2Z).  The  RESULT  is 

pfbf=  Kf'Y-tN0^ 

Now  USE  (2.t)  AND  PUT  £  *  U  .  Thb  RESULT  IS 

. -y  *  r\  2x  i 

p=K<|>  YK(0)JJ\~aoji)  Jfl 

w  0(o) 

Now  WRITE  DOWN  (ZZ)  FOR  -p  ^  DIVIDE  T>  67  p  EXPPESS  Wotf^)  AEftHS 

J  e  J  ^  -irA^O 

OF  (2.6)  aND  (Z.?),  AND  REPLACE  THE  EXPRESSION  THAT  OCCURS, VIZ.  (l-Oo)  , 

y*)  £ 

3Y  Jy'aTs  According  tx>  t3<7)-  The  Result  is 


t*0) 


j=  Y_l  ^  /»  »  %*“•)  . 

Z  -  —  (  0-«oU  )  Ju 

Jb  ”  7-X  a©S  Jo  r  r 


The  INTEGRAL  MAY  BE  INTEGRATED  8|  PARTS  TO  6IVE 

y 


jyyy^ = f  •-«.>* + f~  a°0'-ay>*YJr 

k 

Then  ,  proia  (si),  ( *2),  and  (3l) 

^=\-f 


(S’  •) 


(Si) 


(53) 


WHERE 


t*  i£«  V?> 


(S4) 


The  kinetic  energy  of  the  powder  gas  (41)  may  also  be  expressed  in 
terras  of  S  •'  _ 

(55) 


w,-i£Vt. 


so  that  the  total  kinetic  energy  of  projectile  and  powder  can  be 
expressed  as  9 

MW  =  i  "• v  < '  fO 


<  TOTAL. 


From  (56)  we  see  that  the  total  kinetic  energy  of  projectile  and 
powder  can  be  expressed  as 


AW  =('/2.)me  v 


where  the  effective  mass  mg  is  given  by 


me  =  m  (  \  -t  ~) 


(51) 


(58) 


Wote  that  the  correct  mean  pressure  p  for  use  in  the  equation 
of  state  is  equal  to  the  effective  pressure  ft  that  has  to  be 
assumed  to  act  on  the  effective  mass  Rig  to  impart  to  it  the  total 
kinetic  energy  ('/^ineV.  To  show  this  note  that  Ap£  =  dV/dt 
gives  the  correct  total  kinetic  energy  (on  integration  of  ApeV«At 
Dividing  this  expression  by  Apb  =  m  M/<h  we  findjp  /■»  = 

-  y /j>  as  vas  t0  1:56  shown.  -*b 


) 


Section  VI:  Discussion  of  the  Above  Results 


Letting* 


n  = 


CS9) 


we  have  from  (30^  (3rO>  and  (54) 

5  =  -«cfl*)ndfL  (60) 

6=z(n-n)a.,u (60 

ft- iX’0' •**’>>•  <M) 

In  applications  if  and  6  are  given  andfioand/^  are  required.  We 
first  show  that  the  integral  in  (6f)  can  be  reduced  to  an  ex¬ 
pression  involving  S 


Let 


(63) 


*  n  is  called  the  "polytropic  index"  in  astrophysics, 
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By  expressing  f^^fL  as 

i  r  / 


we  can  express  the  integral  as 


R  '  - 


2(rH-l)£U 


which  can  be  integrated  by  parts  to  give 


_  0-4,1 


1 


I 


R  “  ~ _  _L> 

2Cn-H)<2.p  zCh+Ofl. 


-  J"  a-a-/S)r'(\-a.,/ul)e//i; 


from  which 


2(n+0«..R  =  -(i-«s)MI  +S" 


Solution  for  R  and  division  by  S  gives 

_  _r  _  i  r ,  _  a-*.,)**1 
S'  S  terH^K«L  5 


comparison  of  which  with  (61)  gives 

x  -  -J — r _L  _  Z(^tQ  1 
T  2n+3  La.  €  J 

From  (38)  and  (59)  we  have  also 

*Vfb  ~  (i~^) 


C64) 


From  (64)  and  (65)  it  is  clear  that  the  calculation  of  all  the 
practical  results  of  the  solution  has  now  been  reduced  to  the 
calculation  of  a  from  the  given  values  of  e  and  n  .  We  consider 
this  calculation°and  its  tabulation  in  Section  VII. 


To  set  limits  on  the  parameter  a  we  proceed  as  follows. 

From  (26)  and  (29)  or  from  (38)  we  see  first  that  aQ  must  be  real. 

Now  in  (38)  we  have  ty^Y-l)  =  n*H  and  we  know  that  n  s 
is  a  real  positive  number,  since  y>|.  Thus  if  CL0<O 

we  have  ?./?.  <.  1  ,  an  absurd  result,  since  there  must 

be  a  positive  pressure  difference  to  overcome  the 

inertia  of  the  powder  gas.  If  Q.o>\  ,  may  be 


positive  real  only  if  n  is  an  odd  integer,  any  other  value  of  n 
giving  either  a  complex  or  a  negative  result.  We  may  therefore 
rule  out  the  case  CL«>I  as  non-physical,  since  any  slight  change  of  n 
from  an  odd  integer  value  must  lead  to  a  non-physical  result. 
Thus0ajCU=l  ;  examination  of  (60)  and  (61)  shows 

that  a  =  0  corresponds  to  e  =  0  and  a  =  1  to  e  =  .  For  very 


small  e  we  have,  from  (60)  and  (61),  a  =  e/£2(n+l)].  We  obtain 
the  second  approximation  by  binomial  expansion  of  the  integrand 
in  S  and  of  the  expression  (1 in  powers  of  a  ,  thereby 

obtaining  0‘<lo)‘n“lS  =  ,  +  (l+aJl/j),a-o>  and  then  replacing 

this  aQ  by#  its  first  approximation.  This  second  approximation 
is  cu  rv  &ht3  6 


r\  —  _  & 

2(^  +  0 


D 


LUib 

1 
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insertion  of  which  into  (64)  gives  S'-  3  for  *  O  ,  in  agreement 
with  (f)  and  (3).  Insertion  of  the  first  approximation  for  aQ  into 
(65),  which  may  be  written 

(»-<t©)‘'n  1  —  i  +  (nt»)a0t  - -  ^ 


6  WES 


Jp/f^  -  S'*  k  FaR  € — VO,  IK  AGREEMENT  WITH  (l). 


Section  VII.  The  Methods  of  Calculation  and  Tabulation 


Our  problem  now  is  to  calculate  aQ  as  a  function  of  n  and  e 

from  (60)  and  (61)  and  then  to  devise  tables  by  means  of  which  aQ 

and  1/b  can  easily  bd  found  when  n  and  e  are  given-  i'o  do  so  we 

might  calculate,  for  a  given  value  of  n,  the  quantity  e  for  a 

number  of  evenly  spaced  values  of  a  and  then  inverse  interpolate. 

In  this  way  we  could  construct  a  taBle  of  aQ  versus  e  and  n  for 

eoually  spaced  values  of  both.  In  doing  so  we  should  have  to 

evaluate  the  integral  S,  which  can  be  expressed  as  an  incomplete 

beta-function,  tables  of  which  are  available,  by  Karl  Pearson.  (We 
did,  in  fact,  make  preliminary  computations  in  this  way.) 

In  order  to  explain  most  easily  our  subsequent  calculations 
it  is  desirable  at  this  point  to  discuss  the  tabulation  possibil¬ 
ities.  We  wish  to  construct  tables  with  e  and  n  as  arguments  that 
will  enable  the  user  to  calculate  both  aQ  and  6.  in  these  double- 

entry  tables  it  is  desirable  that  linear  interpolation  be  per¬ 
missible  with  respect  to  both  arguments  s  and  n.  It  seems  point¬ 
less  however,  to  tabulate  both  a  and  b,  since  the  user  of  the 
tables  can  find  either  one  from  tfie  other  by  (64),  undoubtedly 
more  quickly  than  ty  interpolation  in  an  extra  double-entry  table. 

If  the  table  is  of  a  ,  it  is  clear  from  (64)  that  in  calculating  6 
the  user  will  lose  accuracy  by  subtraction,  from  a  table  of  5, 
however,  he  will  be  able  to  calculate  aQ  without  loss  of  accuracy. 

It  is  settled,  therefore,  that  we  must  tabulate  6  or  some  simple 
function  thereof.  The  functions  that  come  to  mind  on  inspection 
of  (64)  are  1/b  and  (2n+3)/5*  •  Since  we  can  probably  insure  the 

permissibility  of  linear  interpolation  with  respect  to  e  only  by 
taking  a  rather  small  interval  for  the  latter  (increasing,  as  we 
see  later,  with  increase  in  e),  our  problem  is  to  find  which  of 
the  three  functions  ^  ,  or(,2n+3)/j  will  give  the  best  inter¬ 

polation  possibilities  for  n.  That  is,  for  fixed  e,  which  of  these 
three  functions  is  most  nearly  a  straight  line  when  plotted  against 
n?  The  above  mentioned  preliminary  computations  show  unmistakably 
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that  only  un+3)A  is  reasonaoly  straight. 

We  showed  above  that  £-3  for6=o  and  that  d9  =  \  corresponds  to 
e—’+oo  ,  so  that  by  (64)  ,  £-*-2X1+3  for  6-**  Oo  .  The  resulting 
simple  table  of  limiting  values 


e 

r 

v% 

(znt 

0 

3 

% 

1+  *"/s 

00 

in +3 

_ 

'/(Zn+5) 

1 

then  shows  that  the  curve  of  VS  versus  n  is  not  expected  to  be 
straight  for  large  values  of  e .  To  see  why  the  curve  of  (*n+S)/f 
versus  n  is  expected  to  be  straighter  than  that  of  6  versus  n 
for  all  e  we  next  examine  the  behavior  of  the  function  (zn+ 3)/s 
in  the  limiting  cases  of  small  n  and  large  fl  .  This  investi¬ 
gation  leads  us  to  the  function  that  we  actually  tabulate. 


To  do  so,  we  first  rewrite  ( 60 )  and  (6l)  in  terms  of  the 
quantity 

(66) 


b  s  a° 


1  -  a0 

From  (60),  (61),  and  (66)  we  then  have 


€=  2(n+l)bjr,£H.b(.-^)j"ju.  (67) 

Since  O  ^  ~  I  y  ve  have  b  ^  0>  b  =  o  corresponding 

to  6=0  and  b“ao  to  e  =«?  *  To  investigate  the  case 

of  small  n  we  write  the  integrand  of  (67)  as  exp{n  InCi-h  bfi-u.*)]  j 
and  expand  this  exponential  in  a  power  series  in  n  .  The  result 
is 


*(n+0b  [i  +  ny'ln  D^b(»-^tjdp-+  ...  J 


Now  as  n  approaches  zero  we  see  that  b  approaches  the  value  e/z 
Any  deviation  of  b  from  the  vaiue6/2,  in  the  logarithmic  term 
will  result  in  a  quadratic  in1)  contribution  to  the  expression 
in  brackets.  Thus,  through  linear  terms  only 


WHERE 


_ _ £ _ x  1  +  nJ"  +  •••• 

2(n+t)b 

=  ( 1 + 1  )'1  In  D + fc  +ltczte)j  ‘lzJ  -1 


(68) 

(69) 
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The  evaluation  or  this  integral  is  given  in  Appendix  c.  r'or  suc- 
f'iciently  small  e  it  has  the  value  e/3. 

From  (64)  and  (.66)  we  have 

z.n+3  „  ,  4.  X  -  (10) 

S  »=>  e 


Comparing  (68)  and  (70),  we  obtain  through  terms  linear  in  n 

-  1  +  2. Jn  +  .... 
y  e 


(71) 


1‘rom  which 


S  =  3  +  rz-  fJ)n  + 


The  initial  slopes  are  thus  2J  for  versus  II  and  Z-*J 

-  e  s  e 

for  0  versus  11  . 


To  investigate  the  case  of  large  n  we  first  note  in  (67) 
that  since  b  $  O  the  integrand  and  thus  the  integral  are  both 
ecual  to  or  greater  than  unity.  Thus  from  (67)  as 

for  a  fixed  value  of  e,  we  must  have  b  -*■  o  .  Since  the  integrand 
then  takes  on  the  indeterminate  form  ,  we  rewrite  it  as  before 

as  exp£n  ln£i  +  b(l-fl*)]J  which  becomes  exp £nb(|-p.*)J 


for  the  limiting  case  b~>  O  •  For  H— *«we  may  also  replace  the 
outside  factor  (rvti)  by  n  ,  thereby  obtaining 


*  =znb  JJ  exp  Cn  b 0  -fi1)]  da 
So  that  b  -  x(6)/fj  1 

WHERE'  K(6)  IS  A  SOLUTION  OF 

c"'V<= 


(13) 

(14) 

(is-) 


Then  from  (70)  and  (74) 

2Jn3  _  ,  +  JL  -  ZfHHl 
s  K(€)  e 

which  gives  for  n  — ►  00 

_l_  _  I _ l_ 

%  "  2K(€)  e 


(76) 

(17) 
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It  is  an,  interesting  check  at  this  point  to  note  thatrt— 
corresponds  to  V— ►  |  ,  i.e.  to  an  expansion  process  governed 

=}’('//>- 1),  which  process  is 

isothermal  if  V ("fc)  is  constant.  The  equation  of  motion  (20) 
can  then  easily  be  integrated  directly,  with  (75)  and  (77)  as 
direct  results. 

On  putting  K  H  S  5  iN  (74)  ,  we  obtain 

e-7r//VH(«-)exp(6'z^  oa.O 

where  ctsk''*  ^  fe* 

and  H(<r)=r  Z(ir",/k)j[  £~3  (H8.Z) 

of  which  there  is  a  very  good  table.*  With  the  aid  of  this 
table  we  may  calculate  the  values  of  e  for  evenly  spaced  values 
of  €~  and,  by  inverse  interpolation,  then  obtain  the  values  of 
and  of  Kstf"2*  for  evenly  spaced  values  of  e. 

Fran  (77)  the  final  slope  of  ?  versus  (1  for  large  n  is 
clearly  zero,  while  fran  (76)  the  final  slope  of  (znfi)/?  versus 

n  is  ( %  -  %)  .  We  may  now  summarize  the  initial  and  final 

behavior  of  the  curves  S'  versus  n  and  versus  fl. 


i  versus  n 

fijpWsn 

Ihitin.  Scope 

z-6J/e 

z% 

Finm.  Slope 

o 

k  -  7e 

Now  J*^£/3  for  sufficiently  small  e  and  J~ &  In  6  for 

sufficiently  large  e,  so  that  the  initial  slope  of  f  versus  n 
varies  with  increasing  e  from  0  to  2  while  the  final  slope  is 
always  zero.  Actual  plots  indeed  show  marked  curvature  for  the 
curves  of  £  versus  n  . 

The  following  table  gives  initial  and  final  slopes  for  the 
curves  (xn+3)/£  versus  H  ,  for  a  few  values  of  e. 


^"Tables  of  Probability  Functions",  Vol.  I,  by  the  WPA,  sponsored 
by  the  National  Bureau  of  Standards  (1941) 


IS 


Slopes  of  the  Curves  Czn+3)/S*  Versus  f) 


6 

Initial  Slope  o< 

Filial  Slope 

Difference 

_ 

0 

0.6667 

0.6667 

0.0000 

2 

0.4929 

0.5602 

0.0673 

•  4 

0.4038 

0.5086 

0.1048 

6 

0.3471 

0.4756 

0.1285 

8 

0.3070 

0.4515 

0.1445 

10 

0.2768 

0.4329' 

0.1561 

0» 

0 . 0000 

0 . 0000 

0.0000 

At  e  -  10  the  difference  of  the  slopes  is  clearly  approaching  a 
maximum  value  which  will  proDably  be  less  than  0.2.  Since  the 
tables  of  the  present  report  end  at  s  =  10,  we  have  not  seen  fit 
to  extend  the  above  table  beyond  that  value,  in  any  case  it 
appears  that  the  curves  of  (Zn+Sj/S*  versus  H  are  very  much 
straighter  than  those  of  $  versus  f)  • 


Further  investigation  then  showed  that  the  quantity Urtt 3)/g 
can  be  expressed  as  follows  in  terms  of  a  very  weakly  varying 
function 


znH.  =  i  +  acnTu- 


(79) 


where  °<  is  the  slope  of  an +3)/r  versus  n  for  n  -o  and  is 

its  slope  for  n=00  .  Thus  ' 

*=  (Qo) 


WHERE 


i T-  yi+T  In  e  +  -  £ 


((,9) 

(81) 
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where  K(€)  is  given  by  (78).  Over  tne  ranges 
Vz  ^n'k  s-  ,  the  quantity  c,  differs  very  little  from  unity, 
ranging  from  0.960  to  1.065.  Such  behavior  is  not  surprising,  for 
a  constant  value  of  C,  will  clearly  reproduce  the  correct  behavior 
of  (Zfl  +3)/$  for  very  small  H  and  for  very  large  11  ;  it  gives 
the  values  \  t(2/fc)  J  for  very  small  n  ,  in  agreement  with  (71)  and 
n  _  zn 

1  +  K(fc)  fc  for  very  large  n  ,  in  agreement  with  (76).  Indeed 
replacement  of  c,  by  unity  in  (79)  gives  values  of  '/£  not  worse 
than  one  part  in  1000  for  any  of  the  values  of  €  and  r)  that  we 
consider.  Such  accuracy  would  probably  meet  the  practical  require¬ 
ments  of  interior  ballistics,  but  it  appears  worth  while  to  obtain 
somewhat  better  accuracy,  since  this  Pidduc£-Kent  solution  has  con¬ 
siderable  mathematical  interest  as  the  limiting  solution  of  a 
classical  problem  in  fluid  mechaUmcs. 

From  (70)  and  (79) 

C,  = 

-3)n 

j  r  i  2.(r\+i)T 

-of  [_  b  e  J 

The  final  C,- table  of  this  report  has  a  range  on  n  from  0.5  to  5, 
corresponding  to  values  of  the  effective  ranging  fron  the  value 
1.2  (approximately  the  value  for  ballistite,  but  uncorrected  x’or 
heat  loss)  up  to  the  value  3,  which  would  allow  for  enormous  heat 
loss . 


C8Z) 

(63) 


We  found  b  by  solution  of  (67),  which,,  for  an  integer  value 
of  n  gives  6  as  a  polynomial  in  b  .  For  II* I  we  had  simply  to 
solve  a  quadratic  equation.  For  n  =  2,  3,  4,  and  5  w-e  tabulated 
£  versus  b  ,  at  equal  intervals,  and  inverse  interpolated  to 
find  b  versus  6  at  equal  intervals.  In  obtaining  a  double¬ 
entry  table  of  C,  versus  £  and  n  we  founa  the  following  argument 
intervals  satisfactory  for  linear  interpolation: 

6  =0(0.2)10)10  AND  n  =  7zf  1,2  34  5 

To  obtain  reasonable  smoothness  in  this  table  we  found  it  necessary, 
because  of  the  subtractions  in  (82)  and  (83),  to  calculate  b  to 
high  accuracy*-,  especially  for  the  smaller  values  of  £  .  This 

circumstance,  conversely,  accounts  for  the  fact  that  very  few 
significant  figures  in  the  C,-  table  will  give  K/%  to  high  accuracy. 


*-  This  necessity  rendered  useless  the  Pearson  tables  of  the  incom¬ 
plete  beta-function. 
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For  fl  =  '/g,  we  ha(i  to  calculate,  according  to  (.67),  values  of 
the  integral 

t'/z  SX'  D+  b<»-/d)3  /lJ^< 

=['/i]l,+ o+b)  (54)** 

In  obtaining  tables  of  b  and  c,  in  this  case  we  found  it  con¬ 
venient  first  to  tabulate  6  versus  b',z'  »  at  equal  intervals 
of  b,/i'  ,  and  then  to  inverse  interpolate  for  the  latter  rather 
than  for  b  . 

We  did  not  have  to  compute  b  or  C»  for  any  other  values  of 
n  ,  except  for  the  purpose  of  checking, which  we  did  very  simply 
as  follows,  -*-he  integral  In s  J<J^i-|.bci-l4t)J0chxin  (67)  satisfies  the 
recursion  formula  f 

,  _  _!__fzfn+i)fi+b)IH  t ) 

fn+i  2n+3L 


We  chose  several  values  of  b  and  for  each  of  these  values  cal¬ 
culated  r„  n-  by  means  of  (84)  and 

successive  applications  of  (85) .  Then  by  (67)  we  easily  computed 
the  value  of  €  for  each  of  these  cases  ana,  from  (70),  the 
accurate  value  of  */S  .  Tc  check  the  adequacy  of  the  unit  rt-interval 
of  the  .  C,- table  we  first  interpolated  linearly  in  the  latter  (or 
used  the  gxaphs  or  contour  map  given  later).  We  then  interpolated 
in  the.,  tables  of  ot  and  £  to  find  the  values  of  the  latter 
quantities  corresponding  to  each  £  and  then  inserted  the  values 
of  n, ana  c,  into  (79)  to  find  */£  .  Comparisons 

with  the  accurate  values  of  '/S*  showed  that  the  error  was  never 
reater  than  that  arising  from  the  use  of  tabular  values  of  n  . 
See  Section  VIII  for  a  discussion  of  errors  of  the  tables.) 

For  the  table  of  we  computed  anchor  values  directly  from 

the  formula 

“  I[0tf) /z  Inlj+fc*  Je(z7ej  -2j  (st>) 

**  For  the  derivation  see  Appendix  C. 
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deducible  directly  from  (69)  and  (80).  These  anchor  values  were 
6  =  0(OM)2, 6(o. SL)s'.o(o.r)  10.0.  Direct  interpolation  then 

gave  values  of  for  the  final  tabular  arguments-:  6-0(0.05)  1.0 

(0.1)  4*0.  (0.2)  10.0. 

For  the  table  of  /&(£>  we  first  used  (78),  wi  th  the  aid  of  the 
WPA  "Tables  of  Probability  Functions",  to  make  a  table  of  fe  for 
ff-  *  0.18  (0.02)  1.30.  Inverse  interpolation  then  gave  tables  of 
and  thus  of  for  the  same  anchor  values  of  £ 

as  for  cl  .  Division  gave  a  table  of  ^  for  those  same  anchor 

values  and  direct  interpolation  then  gave  the  table  of  7^  for  the 
final  tabular  arguments:  €.=  0(0.1)  2.6  (0.2)  5*0  (0.25)  10.0. 


Section  VIII.  Discussion  of  the  Tables  and  their  Accuracy 

Appendix  A  gives  the-  final  working-  tables  for 
and  C,  (6,n)  .  Appendix  h  gives  a  set  of  graphs  of  c,  versus 

£  FdR.  r?  s  '/z. ,  i>  t,  +,  5*  ,  a  contour  map  of  C,  as  a  function 
of  6  and  rj  and  a  table  of  values*  read  from  the  graphs  and 
used  in  constructing  the  contour  map. 

The  values  of  cL  and  ^3  are  given  to  four  decimal  places  and 
are  correct  in  all  cases  to  four  decimal  places,  the  original  anchor 
values  of  <k  and  ©^having  been  calculated  to  seven  decimal  places. 
The  values  of  C,  are  given  to  three  decimal  places,  rounded  down 
from  four  decimal  places.  The  possible  error  in  a  tabular  value  of 
°c  or  of  ^  is  thus  not  greater  than  0.5  in  the  fourth  decimal,  i.e 
0.00005.  Examination  of  the  differences  shows  that  linear  inter¬ 
polation  is  permissible  throughout  either  table.  With  a  maximum 
possible  erro^  of  absolute  value  0.00005  in  either  of  two  neighbor¬ 
ing  entries,  it  is  clear  that  the  permissible  linear  interpolation 
cannot  give  an  error  of  magnitude  greater  than  0.00005  either  for  oc 
or  for  .  Experience  with  the  c,  table  (which  is  not  subdivided 
quite  finely  enough  to  make  linear  interpolation  always  rigorously 
permissible)  and  with  the  C,  graphs  and  contour  map  shows  that  the 
possible  error  is  one  unit  in  the  third  decimal  place,  i.e.  0.001. 

We  now  ask  the  question:  how  much  error  can  be  produced  in  '/S' 
by  these  possible  errors  o.oooof.  =  o.oooos  e,ool  7 

We  have  r  J  ' 

ICti  =  i  +  ^nr  ,n9) 

$  L  i  +  c.n  J 


*  This  table  is  included  to  facilitate  possible  reproduction  of  this 
report  in  other  forms. 
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A  word  of  warning  is  appropriate  here.  In  evaluating  the  fraction 
O+c, it  would  be  entirely  incorrect  to  keep  only 
three  decimals  in  the  numerator  and  three  in  the  denominator,  simply 
because  c,  is  given  only  to  three  decimals.  One  can  best  see  this 
point  by  considering  tile  case  e=zln=»  for  which  1.1364 
and  C,  =■  1.061.  We  compute  the  fraction  {1+ +  as 

1.070  for  C,  *1.061  and  1;068  for  C,  *  1.  Thus  rounding  c,  to  zero 
decimal  places  leaves  unaffected  the  second  decimal  place  of  the 
fraction,  so  that  we  are  justified  in  keeping  two  more  decimal 
places  in  the  result  than  were  used  for  C,  .  Instead  of  using 
elaborate  rules  for  the  proper  number  of  significant  figures, 
however,  one  will  probably  find  it  simpler  to  use  more  places  than 
are  necessary  or  valid  in  calculating  (79)  and  simply  round  down 
at  the  end  to  the  proper  number  of  figures  for */£  that  are  in¬ 
dicated  by  the  considerations  about  to  be  given. 


Error  in  v%  from  Error  in  cC  Alone 

Let  A*  e/s)  be  the  error  in  Vt  arising  from  insertion  into  (79) 
of  a  value  of  ©<  in  error  by  .  From  (79)  we  have 


Placing  C,  -/  and  (its  value  for  6=  10)  we  find  that  the 

coefficient  has  its  largest  value  for  5  (for  the  range  of  values 
of  n  that  occur  in  the  tables).  This  value  being  0.565,  we  find 
for  thfe  largest  possible  error  in  '/f  that  can  arise  from  the  possible 
error  0.00005  in  aC  t  the  value 

^  0/S)_  Q.OOOOS. 

or  three  units  in  the  fifth  decimal  place  of  '/£  . 


Error  in  Vs  from  Error  in  ^  Alone 

Let  &f('/S)  be  the  error  in  '/$  arising  from  insertion  into  (79) 
of  a  value  of  4  in  error  by  a/?  .  Placing  C,=  1  in  (79)  we  then 
have  ' 

A 


‘f)- 


OCX) 


(n-h)fzn+3) 


V* 


(88) 


The  largest  value  of  n  (i-H))",(zn+3)  for  our  range  of  values  of 

n  occurs  at  n=5  and  has  the  value  0.32.  The  largest  value  of  ©c 
is  V3  for  f-  =  0.  For  a  possible  error  ~  0.00005  we 

thus  find  ^  r 

=  fo.3Z)(|)(o.OOOO$-)  -o.oooo 
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or  one  unit  in  the  fifth  decimal  place  of  ‘/r. 


Error  in  Vi  from  Error  in  C,  Alone 


Let  ( '/£)  be  the  ffror  in  Vs  arising  from  insertion  into 
(79)  of  a  value  of  C,  in  error  by  &C(  .  From  (79)  we  have 


AC|(  VS) 


on*  ac, 

0+n)z,Czn+3) 


(09) 


on  differentiating  (79)  and  putting  •  The  largest  value  of 

ot(fg-\) is  about  0.16,  occurring  for  €.  =  10,  and  the  largest  value 
for  the  function  i  +  n)*'^fZT^^3)_,  is  about  0.064,  occuring  at  t\Cc 

2.3.  On  insertion  of  these  values  and  the  maximum  possible  value 
0.001,  we  find 

oiO 

Ac  -  0, 0000 1 

MAX  ’  WV| 

or  one  unit  in  the  fifth  decimal  place  of  ‘/f. 

Thus  in  the  practical  use  of  Eq.  (79)  and  the  tables  for 
and  C,  we  may  expect  a  maximum  possible  error  units 

in  the  fifth  decimal  place  of  l/f  •  Now  the  smallest  Vs  for 
our  ranges  of  £  and  h  ,  occurring  for  €=10  and  11  =  5,  is  about 
0.23,  so  that  the  maximum  possible  error  is  about  5  parts  in 
23,000  or  1  part  in  46QO.  Roughly,  then,  the  maximum  possible 
error  is  about  1  part  in  5000  in  calculating  •/$  by  (79)  and  the 
tables.  It  would  be  quite  feasible  to  reduce  this  possible  error 
to  a  much  lower  value  by  using  finer  intervals  and  more  decimals 
in  the  tables  of  oC  and  ^  ana  still  retain  the  desirable  feature 
of  permissible  linear  interpolation.  With  the  use  of  such  bulkier 
tables  of  ©*  and  one  could  reduce  the  possible  error  to  that 
arising  from  error  in  C,  alone,  or  about  1  part  in  25,000.  We 
have  not  felt  it  worth  while  to  do  so,  however. 


Possible  Error  in  Up 


From  (64)  we  have 

_L  _  2.n+3  z£rvH) 

ao  ~  £  T  £ 

so  that 

Aa0  -  ~<20z£zn+3M  (  j) 


From  (90)  and  (91) 


Aa0 

do 


zn+3 


^n+3  z(n-H) 
%  *~e 
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(90) 

(9\) 

(9Z) 


We  have  seen  that  by  far  the  largest  part  of  the  possible  error 
in  y$  arises  from  the  possible  error  in  <*  and  that  the  largest 
possible  error  in  '/f  aue  to  error  in  occurs  at  h  =  5.  To  obtain 
an  upper  limit  to\*<Lo/Q.0\  ,,e  therefore  insert  into  (92)  the 

value  n  =  and,  obviously,  the  largest  value  of  6  ,  viz.  10,  that 
occurs  in  our  tables.  Using  the  value  (*0+3)/?  as.  3.034  for 
6-10  and  H-  5  and  the  value  I  A  ( i/f ) I v  0.00005  we  find 


The  maximum  possible  error  in  4,»  is  thus  only  1.5  parts  in 
10,000  or  1  part  in  6700.  It  is  clear  that  the  maximum  possible 
absolute  error  also  corresponds  to  the  same  values  of  €  and  tj 
for  which  0.2363,  so  that  laapl,^  =  0.000035. 

Possible  Error  in-^/jb 


The  Pidduck-Kent  value  of  the  ratio  of  breech  pressure  to 
Dase  pressure  is  given  by 


Thus 


-  n-i 


^(F0/fb)  »  (nti)O-O 


us) 


and 


A<P«,/;Pb)  (n+Qatto 
(■Po/PJ  "  '•*- 


(93) 


On  inserting  the  values  A  a0  -  0.000035,  ao  =  0.2362,  and  n=  5. 
ve  find 


wax  I  A(j>./rk) 

I  (p./h) 


=  o.ooo  zn 


The  maximum  possible  error  that  can  arise  in  calculating  the 
Pidauck-rKent  value  for  this  pressure  ratio  from  our  tables  is 
thus  2.7  parts  in  10,000  or  1  part  in  3700. 

Possible  Error  in^/Jb  OR  /W-p 

The  Pidduck-Kent  value  for  the  ratio  of  space  mean  pressure 
to  Dase  pressure,  or  for  the  ratio  of  the  total  kinetic  energy  of 
the  powder  plus  projectile  to  that  of  the  projectile,  is  given  by 


The  greatest  possible  error  in  either  ratio  is  thus  £  times  the 
greatest  possible  error  in  '/£  or  0.00005  6  .  For  €.=■  10  and  n-5 
the  fractional  error  thus  becomes 


0.0005 _ 

1  +  (0.~23J  10 


0.0002 


The  greatest  possible  error  that  can  arise  in  calculating  the 
Pidduck-Kent  value  for  this  ratio  from  our  tables  is  thus  1  part 
in  5000. 

v  < John  P.  Vinti 


Sidney  Kravitz 
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APPENDIX  A 


TABLES  POR  THE  PIDDUCK-KENT  SOLUTION 

S  zn+i\  l  i+ 


where  n  =  l/(  Jr  -  1) 
1‘ollow  below. 


(12±£L\\ 

\  (+c,  n  JJ 


and  <*, 


,  and  C,  are  given  in  the  tables  which 
TAdlE  I 


In  the  following  table  ©t  is  given  as  a  function  of  €  for  the 
following  values:  £  =  0( .05)1( .1)4( .2)10.  It  is  found_£rom  the 
following  formula:  ^  =  ±fyTTY-[ln  0  +  6 +6JT7J )}  -3 

The  first  forward  differences  of  the  function  are  given  in 


the  third  column  marked  A].  Linear  interpolation  is  permissible. 


€ 

A, 

€. 

A,  j 

£ 

©c 

A, 

.00 

.  6667 

-66 

1.0 

.5621 

-81 

4.0 

.4038 

-67 

.05 

.6601 

-64 

1.1 

.5540 

-78 

4.2 

.3971 

-65 

.10 

.6537 

-62 

1.2 

.5462 

-75 

4*4 

.3906 

-62 

.15 

.6475 

-61 

1.3 

.5387 

-73 

4.6 

.3844 

-59 

.20 

•  6414 

-59 

.1,4  _ 

.5314 

-70 

4.8 

.3785- 

-57 

.25 

.6355+ 

-57 

1.5 

.  5244 

5.0 

.3728 

-$"6“ 

.30 

.6298 

-56 

1.6 

.5177 

-65 

5.2 

.3672 

-53 

•  35 

. 6242 . 

-55 

1.7 

.5112 

-63 

5.4 

.3619 

-51 

.40 

.6187 

-53 

1.8 

.5049 

-61 

5 . 6 

.3568 

-49 

.45 

.6134 

-52 

1.9 

.4988 

-59 

5*8 

.3519 

-43 

.50 

~imr- 

-51 

2.0 

.4929 

-57 

6.0 

.3471 

-46 

.55 

.6031 

-50 

2.1 

.4872 

-56  | 

6.2 

.3425+ 

“44 

.  60 

.5981 

-43 

2 . 2 

.4816 

-53 

6  e  4 

•  3331 

-43 

.65 

.5933 

-48 

2.3 

•  4763 

-53 

6 . 6 

.3338 

-42 

.70 

. 5885+ 

-46 _ 

2*4 

.4710 

-50 

6.8 

.3296 

-41 

.75 

.5839 

-46 

2.5 

.  466O 

-50 

"770 

.3255+ 

.80 

.5793 

-44 

2.6 

.4610 

-48 

7.2 

.3216 

7-7  iil 

.85 

-44 

2.7 

.4562 

-46 

7.4 

.3178 

-37 

.90 

.5705+ 

-42 

2.8 

.4516 

-46 

7.6 

.3141 

-36 

•  95 

.  5663 

-42 

2.9 

.4470 

-44 

7.8 

.3105+ 

-35 

3.0 

•  4426 

“4^ 

“570 

. 3U7U — 

3.1 

•  4383 

-42 

8.2 

.3036 

-33 

3.2 

.4341 

-41 

8.4 

.3003 

-32 

3.3 

-40 

8.6 

.2971 

-31 

3.4 

.4260 

-39 

8.8 

.  2  9  A0 

-31... 

3.5 

.4221 

-39 

9.0 

.  2909 

-30 

3.6 

.4182 

-37 

9.2 

.  2879 

-29 

3.7 

*  4145  + 

-36 

9.4 

.2850+ 

-28 

3.8 

.4109 

-36 

9.6 

.2822 

-27 

3.9 

.4073 

-35 

9.8 

.2795- 

-27 

10.0 

.2768 
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TA8LE  II 


In  the  following  table  jS  is  given  as  a  function  of  6  for 
the  following  values:  €.  =  0(,1)2.6(.2)5(*25)10. 

I  rl  _  2.1  k  f  1  -KM 

P  ~  S?  L  K  "  where  k  is  the  solution  of  €  —  2  K  C.  J0 


The  first  forward  differences  of  the  function  are  given  in 
the  column  marked  Ai.  Linear  interpolation  is  permissible. 


e 

- 

A  * 

6 

/* 

A. 

€  " 

A  _ 

A, 

0.0 

1.0000 

67 

2.5 

1.1686 

63 

5.00 

1.3160 

136 

0.1 

1.0067 

69 

2.6 

1.1749 

125 

5.25 

1.3296 

135 

0.2 

1.0136 

69 

— 

— 

— 

5.50 

1.3431 

134 

0.3 

1.0205- 

69 

2.8 

1.1874 

124 

5.75 

1.3565+ 

132 

0.4 

1.0274 

70 

0.5 

1.0344 

69 

3.0 

1.1998 

122 

6.00 

1.3697 

131 

0.6 

1.0413 

70 

3.2 

1.2120 

121 

6.25 

1.3828 

129 

0.7 

1.0483 

70 

3-4 

1.2241 

119 

6.50 

1.3957 

128 

0.8 

1.0553 

69 

3.6 

1.2360 

118 

6.75 

1.4085- 

126 

0.9 

1.0622 

69 

3.8 

1.2478 

117 

1.0 

1.0691 

69 

4.0 

1.2595- 

115 

7.00 

1.4211 

126 

1.1 

1.0760 

69 

4.2 

1.2710 

114 

7.25 

1.4337 

124 

1.2 

1.0829 

68 

4.4 

1.2824 

113 

7.5Q 

1.4461 

123 

1.3 

1.0897 

68 

4.  6 

1.2937 

112 

7.75 

1.4584 

121 

1.4 

1.0965+ 

68 

4.8 

1.3049 

111 

1.5 

1.1033 

67 

8.00 

1.4705+ 

121 

1.6 

1 . 1100 

67 

8.25 

1.4826 

120 

1.7 

1.1167 

66 

8.50 

1 . 4946 

118 

1.8 

1.1233 

66 

8.75 

I.5O64 

118 

1.9 

1.1299 

65 

2.0 

1.1364 

66 

9.00 

1.5182 

116 

2.1 

1.1430 

64 

9.25 

1.5298 

116 

2.2 

1.1494 

65 

9.50 

1.5414 

115 

2.3 

1.1559 

64 

9.75 

1.5529 

113 

2.4 

1.1623 

63 

10.00 

1.5642 
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TABLE  III 


In  the  following  table  is  given  as  a  function  of  e  and  n 

for  all  combinations  of  the  following  values:  e  =  0(.2)l(l)l0; 
n  =  1/2, 1(1)5.  . 


-  n  — 1 

e 

1/2 

1 

2 

3 

4 

5 

0.0 

1.000 

1.000 

1.000 

1.000 

1.000 

.2 

1.016 

1.016 

1.016 

1.016 

1.016 

1.016 

.  029 

1.029 

l-ffl 

1.046 

1.029 

1.029 

1.030 

:I 

- 

L.Q38 

L.045 

1.039 

1.046 

1.039 

1.047 

1.039 

1.047 

1.039 

1.047 

i 

1.051 

1.051 

1 . 052 

1.052 

1.053 

1.053 

1.059 

1.061 

1.063 

1.064 

1.065 

1.065 

3 

1.053 

1.057 

1.061 

1.063 

I.O64 

1.065 

4 

•  1.042 

1.047 

1.053 

1.055 

1.057 

1.058 

5 

1.029 

1.036 

1.042 

I.O46 

1.048 

1.049 

6  . 

1.015 

1.023 

1.031 

1.035 

1.037 

1.039 

7 

1.000 

1.010 

1.019 

1.024 

1.027 

1.029 

8 

.986 

.997 

1.008 

1.013 

1.016 

1.018 

9 

.973 

.985 

.997 

1.002 

1.006 

1.008 

10 

.960 

.973 

.986 

.992 

.996 

.993 

The  above  table  is  represented  graphically  in  Graph  I  of 
Appendix  B  which  follows.  Appendix  B  also  contains  a  contour 
map  of  c,(e,n)  (given  in  Graph  II),  and  the  tabular  values 

used  to  plot  Graph  II  (given  in  Table  IV) . 
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APPENDIX  B 

GRAPHICAL  REPRESENTATION  OF  THE  TABLE  FOR 
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TABLE  IV 


In  Graph  II  on  the  preceding  page  n  is  plotted  as  a  function 
or  e  for  constant  values  of  c^.  In  the  following  table  the  values 

used  to  plot  these  curves  are  tabulated  with  e  given  as  a  function 
of  c^  and  n. 


_ : _ _ _ n  _ I 

C1 

1/2 

1 

2 

3 

4 

5 

1 .  GOO 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.010 

0.11 

.0.11 

0.11 

0.11 

0.11 

0.11 

1.020 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

1 . 030 

0.43 

0.42 

0.42 

0.41 

0.40 

0.40 

1.049 

0.65 

0.65 

0 . 64 

0.64 

0.61 

0.60 

1.050 

0.97 

0.97 

0.94 

0.92 

0.91 

0.90. 

1 . 0.60 

— 

1.69 

1.50 

1.44 

1.40 

1.35 

1 . 063 

— 

— 

1.98 

1.75 

1.63 

1.60 

1.065 

— 

— 

— 

— 

2.00 

1.92 

1.065 

_ _  _ 

_ _ 

_  __ _ _ 

2.73 

2.98 

1.063 

— 

— 

2.45 

2.95 

3.18 

3.35 

1.060 

— 

2.50 

3-12 

3-43 

3.6? 

3.75 

1.050 

3.29 

3.73 

4.23 

4-55 

4-75 

4.87 

1.040 

4*15 

4 . 60 

5.20 

5-52 

5.73 

5.89 

1.030 

■  4-92 

5-44 

6.06 

6.45 

6.70 

6.35 

1.020 

5-65 

6 . 23 

6.90 

7 . 37 

1  7.65 

7.85 

1.010 

6.34 

6.98 

7.80 

8.30 

|  *‘60 

8.83 

1.000 

7.00 

7.75 

8.70 

i  • 

\  9  *  5  ^  J 

9.80 

.990 

7.72 

8.57 

9.65 

— 

— 

.980 

8.45 

9.40 

— 

— 

— 

— 

.970 

9.20 

— 

— 

— 

— 

— 

APPENDIX  C 

EVALUATION  OF  CERTAIN  INTEGRALS 


(The  auxiliary  variables  used  in  this  appendix  have  no 
connection  with  "the  notation  of  the  rest  of  the  report.) 

1.  Evaluation  of  JT  S  f0  Q 

On  putting  I  +  ^  %.)  ^  AtfO  (&/z)  SX|^ 

vefind  J-f|  j'Zf*'^(K'-XZUX. 

=  ( I)''"//'  +  'Ack-*i^* 

-  (  j)VTKA>££+  X,  A-fK^-X,1) -2x0 

,  T  T  Ktx.  _  KLf  Xi1*  +  iKX,  _  ^  :  -r 

Mew  K2--/,1 -t  AUD  X-X,'  "kv-x7 

thus  J*  *  (i  t|r  in  Vect^TJ  - 1 

Z>  Evaluation  of  X'/2S-^  C*  +  bO-xt)J 

^  iXi  -o+b)^/:c«  -^j^ 

Ler  X  C  M* b)  s  ^  THAT  -&■  RA^GS  FROM  ^  TO 

-  .S'»~'tThy*  =  +an’'  b 

_  —  V  .  -I/-. 

I  HEN 

Nd*J 

Thus 


X7l  =0+b)  b'1'*-//"'  <**'+*» 


f*m  C0&le-cH&  =('/?.)(&m  +  SlhA»l 


I^S  (fcjfT+D+uV'^W 


••to 


3.  Qiv€N  Ih  *  jV  Dl"  b(l-x*-)]ndx  ^  TO  DERIVE  A  RECURSION  FORMULA 


FOR  T, 


n+l  /.l  ^ 

I„  ,-J.  Cl  +  b(i-x*-)jn C Hr b( I -x*-XJ dx 


In+1-(l+b)ln  =  -b  j  £L+b(l-x*jJ  n  x2dx 

=  1/2  (  x£l+b(l-xB)J  n  d[l+b(l-xs)] 

Xro 

=  2(n+l)  ^l+b(1^x2)]n  1\o  ~  J^+l)f  [1+b^1~x2]]  n  dx 

*0 

_  1 _  In+T. 

2(n+i;  2(n+l) 

*“»  ■'  fill)  In+1  ‘  tUb)In  +  irW 

“>»  w  ■  sfc-[>(i»i)(i*wiB  +  0 

We  therefore  have 

b/2  =  (V4)[3(l+b)I1/2  .  l] 

X5/2  =  (l/6)[5(l+b)I3/2  +  l] 

h/2  =  (l/8)[7(l+WI5/2  +  l] 

i9/2  =  (i/io) [9(1+ b)i7/2  t  r] 
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TABLE  OF  svMun^ 


Page  where 
symbol  first 
appears 


uniform  cross-sectional  area 
a  parameter 

a  parameter  characteristic  of  the  Pidduck-Kent  solution 
a  constant 

V(iv 

initial  mass  of  the  powder 
a  function  of  e  and  n 
specific  heat  at  constant  volume 
total  internal  energy 
a  function  of  x  alone 
abbreviation  for  f ' (x) 
abbreviation  for  f(x^) 
probability  integral 

=  yo'c+  b<i-xi>jndx 


-Jo'  /nCi  +  | 

a  solution  of 


( t -/£*)]  Jll 

*  -  *KeK/.’e 


mass  of  the  projectile 
effective  mass 

the  polytropic  index,  1/(8  -l) 
pressure  of  the  gas 

pressure  at  the  base  of  the  projectile 
effective  pressure 
breech  pressure 


10 

3 

10 

9 

16 

3 

19 

( o 

II 

7 

8 

10 
18 

21 

16 

17 

3 

13 

13 

5" 

3 
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T> 

Q 

R 

S 

T 

t 

V 

\ 

V, 


space  mean  pressure  behind  the  projectile 

>  r  J _ 2(n-H)  ~| 

~  Lb  £  J 

the  gas  constant  per  unit  mass 

absolute  temperature  of  the  powder  gas 
time 

velocity  of  the  projectile 

kinetic  energy  of  the  powder  plus  powder  gas 
kinetic  energy  of  the  projectile 


Vi  total  total  kinetic  energy  of  projectile  and  powder 


w 

vc 

X 


the  variable  { 


P-1) 


z 

oC 

f 


the  variable  W  <|» 

initial  distance  of  any  gas  particle  from  the  breech 

initial  distance  from  the  breech  to  the  base  of  the 
projectile 

distance  from  the  breech  at  time  t  of  a  particle  whose 
breech  distance  was  initially  x. 

reduced  breech  distance  of  a  particle 

initial  slope  of  (2n  +  3)/b  versus  n. 

ratio  (final  slope/ initial  slope)  of  (2n  +  3)/b 
versus  n 


effective  ratio  of  specific  heats 
y*  ratio  of  specific  heats  of  the  gas 

S'  =  S/R 

6  charge-projectile  mass  ratio 

specific  covolume  of  the  powder  gas 
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!*-  dummy  variable  iO 

f  density  of  the  gas  5 

average  density  U 

fo  initial  density  of  the  gas  in  the  special  solution  5" 

G-=  kV2 

16 

functions  of  time  alone  6 
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